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ABSTRACT 

We analytically investigate gravitational radiation induced by a test particle 
falling into an extreme Kerr black bole. Assuming tbe radiation is dominated 
by the infinite sequence of quasi-normal modes which has the limiting frequency 
m/(2M), where m is an azimuthal eigenvalue and M is the mass of the black 
hole, we find the radiated energy diverges logarithmically in time. Then we 
evaluate the back reaction to the black hole by appealing to the energy and 
angular momentum conservation laws. We find the radiation has a tendency to 
increase the ratio of the angular momentum to mass of the black hole, which 
is completely different from non-extreme case, while the contribution of the test 
particle is to decrease it. 
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1. Introduction 


Coalescence of binaxy neutron stars or black holes has been considered as one 
of the strongest sources of the gravitational radiation. As for the total energy ra- 
diated in the process of the coalescence, several analyses have been done. Gilden 
and Shapiro [ll calculated the head-on collision of two neutron stars by Newto- 
nian hydrodynamics code and obtained the typical efficiency (A E/M) of 0.1% 
for colliding 1AM Q neutron stars. As for fully general relativistic calculations, 
Smarr l2l performed simulations of the head-on collisions of two black holes and 
obtained the efficiency of ^ 0.1%. Stark and Piran 1,1 calculated the formation 
process of rotating black holes and obtained the efficiency of ^ 0.1%. One of 
the important results in these calculations is the wave pattern is characterized 
by a quasi-normal mode (QNM) of the black hole except for the very early stage. 
Especially the wave pattern by Stark and Piran ^'is very similar to that from 
a perturbation calculation of a rotating dust ring of mass /z falling into a Kerr 
black hole of mass M (^> /z)! 41 The collapsing dust ring mimics the flattening 
in the equatorial plane of the collapsing polytrope. One also finds very good 
agreement in the amplitudes if one scales /z up to the reduced mass of the system 
which is no longer infinitesimal. These experiences suggest us that the pertur- 
bation calculations are a very good guide to fully nonlinear general relativistic 
calculations!* 1 

One may think that the above fully nonlinear simulations show the weakness 
of the sources of gravitational waves. However all these calculations are restricted 
to axially symmetric systems. Perturbation calculations for the efficiency of the 
emission of gravitational waves from a test particle falling into a Schwarzschild or 
Kerr black hole for various cases l4) show that in general the efficiency for axially 
symmetric systems is much smaller than that for non-axially symmetric systems 
due to the phase cancellation effects!* 1 The perturbation calculations for Kerr 
cases [4,T '* 1 show that the wave pattern is also dominated by a QNM of the black 
hole which depends on the azimuthal eigenvalue m. The energy from a particle 
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with orbital angular momentum falling into a Kerr black hole increases with the 
increase in the Kerr parameter a. This is due to the fact that the imaginary part 
of the fundamental QNM for m = 2 goes to zero in the limit of a = M where M 
is the mass of the black hole. The small imaginary part means the large damping 
time, hence much energy. However what will happen for an extreme Kerr black 
hole (o = M) case was not clear in the previous studies since they were only 
for a < 0.99 M. In this paper we shall study analytically the gravitational waves 
emitted from a particle falling into an extreme Kerr black hole as a first step to 
know the maximum efficiency of the gravitational radiation. 

One may think that the formation of an extreme Kerr black hole is of only 
academic interest. But it is not so. Observationally there are 10 binary pulsars 
at present. Eight of them have small eccentricity and small mass companion 

0.4Af Q ). Three of them have large eccentricity and large mass companion 
(~ 1.4M q ). Especially PSR1913+16 is believed to consist of two neutron stars 
of mass 1.445Af Q and 1.384M 0 . [#I Due to the emission of gravitational waves, two 
neutron stars in PSR1913+16 will coalesce in ~ 10® yrs. The recently discovered 
binary millisecond pulsar PSR0021-72A [lol will coalesce in much shorter time ~ 
10® yrs. Under the assumption of a steady state, it is estimated that coalescence 
of binary neutron stars will occur ~ 10 events/year up to a distance of 100 Mpc. 
Therefore they can be relevant sources of gravitational waves. Now the total 
mass of such a binary like PSR1913+16 is much larger than the maximum mass 
of the neutron star. Therefore a natural final destiny of such a coalescence is the 
formation of a black hole. Moreover the total angular momentum of the system 
when the coalescence begins is larger than that for the extreme Kerr black hole. 
So the final black hole can be an extreme Kerr black hole. 

This paper is organized as follows. In §2, we review the property of QNMs of 
the extreme Kerr black hole. In §3, the wave form of the gravitational radiation 
induced by a test particle falling into a black hole is analytically estimated. 
The effect of back reaction is also considered. §5 is devoted to conclusions and 
astrophysical implications of the results. 
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2. QNMs of the extreme Kerr black hole 

There exist infinite number of QNMs which accumulate onto the critical 
frequency u> c = am/(2Mr + ) = m/(2M) for a = M. This occurs for modes with 
\m\ = l and it has been shown both analytically by Detweiler lu 'and numerically 
by Leaver. Following Detweiler, we first briefly review the derivation of this 
fact, rewrite the result in the form more convenient for us, and then discuss the 
properties of the QNMs. 

The QNMs correspond to complex zeros of Z- n ( u>), the amplitude of ingoing 
radiation at infinity of the mode which is purely ingoing at the horizon, called the 
”in” mode. Its form at a ~ M has been discussed in Appendix A of Teukolsky 
and Press! 15 ' Taking the limit a — ► M, one finds that 

z m = a(-i,ii)e i '(' +3 '' 2 )(2ii)- !, -‘(-8a.T) ,+1 / 2 - i (“ +s ) fre" - , 

( 2 . 1 ) 

where 


a? = u)M , 
m 

T ~ u 2 ’ 


a(-£,a>) = 


r(-2iS)r(l - 2 iS) 


T(s + 1/2 — x( 2u> + S))T( — s + 1/2 — i(2Q + £)) 


r = 


a(6,u>) 


a(-6,u) 


cosh 7r(2a> — 6 ) 
cosh7r(2u> + £) ’ 


( 2 . 2 ) 


8 = arg 

i 

S 2 =m 2 -(s + i) 2 - J A, 

with s being the spin of the wave (a = ±2 for gravitational waves) and ,A, the 
eigenvalue of the spin- weighted spheroidal harmonic which depends on d>, l and 
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a(-£,d>)| 



m (for notational simplicity, we omit the angular indices ( l,m ) unless there is a 
chance of confusion). For our purpose, we shall take s = —2 and suppress the 
spin indices in the following. In addition, it is known that S is real for the case 
of our interest (i.e., for all modes with |s| = 2 and |m| = / > 2, where m > 0 for 
Re(u>) > 0 and m < 0 for Re(u;) < 0)! l>,,4) Then it can be taken as a positive 
number without loss of generality. 


For later convenience, we translate the above formula for Z. into the am- 
plitude A in of the corresponding ”in” mode solution for the equation derived 
by Sasaki and Nakamura| 1&1 which has a nice property of being regular at both 
infinity and horizon. The derivation of the relation between Z in and A. is given 
in the Appendix. One finds 


A 


m 




(2.3) 


Since the infinite sequence of QNMs, which is given by solving for zeros of 
the square bracket term in Eq.(2.1), appears for u> very close to u> c = m/2 (».e., 
for r ~ 0), one may replace <2> by u> c in the rest of the expression (2.1) or (2.2). 
Letting — u>r = z, one obtains after a straightforward algebra that 


A. («;6.) = e » coshlr j 2 ^ + *) 
,nV c ' iu) c smh 2ttS 


re i0 ° - z 2iS 


where 


= arg[a(-^,d> e )] - (2a> c + S) In 8 , 

0 Q = 0 — 261n8. 

Then it is readily seen that the zeros of A. are given by 


(2.4) 


(2.5) 


flo— 2 t»t 


Z = z n = C 26 


Wc 


( 2 . 6 ) 


for any (large) positive integers of n, where the phase <p c is given by 

*, c = -ii nr 

cosh7r(2u’ c 5) 
cosh7r(2w c — ^)J 



(2.7) 
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In terms of u>, these QNM frequencies axe 


1 flo — 2ny 


u n =u c + T n =u c - — e 22 


To see the properties of <2> n in more details, let us rewrite the above result by 
noting the expression (2.7) for the phase y? c and the fact that <2> e = m/2. We 
then obtain 


"n = y - 

2 771 


2 0o— 2nir 2 0o — 2nr 

~ e 26 cos IV’cl -*T— T e 25 sin|p c l» 

[ml 


where 


I,. | _ J_ i_ [ cosh7r(|m| +f) ] 

^ c ' 25 1 [cosh7r(|m|-^)J 

Thus the imaginary part is always negative, irrespective of the sign of m. Note 
also that the QNMs of Eq.(2.9) have the symmetry u> n</ _ m = In any 

case, this implies that an extreme Kerr black hole is stable against the excitation 
of a finite number of the QNMs belonging to the infinite sequence given above. 
However, since the imaginary part tends to zero as n approaches infinity, it is 
not obvious that the excitation of an infinite number of the QNMs would not 
lead to some kind of instability. In fact, we shall see in the next section that the 
radiation induced by a test particle diverges logarithmically in general. 


3. Gravitational radiation induced by a test particle 

In this section, we investigate analytically the wave form of gravitational 
radiation induced by a test particle of energy e falling into an extreme Kerr 
black hole of mass M. We assume e <C M. 

The basic equation we consider is the radial equation derived in Ref.15 which 
takes the form, 

where the explicit forms of T and W, and the relation of the source term S to 
that of the original Teukolsky equation are given in Ref.15. Here we only mention 
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that T and U vanish sufficiently rapidly for both r* = ±oo, and so does 5 for 
any trajectory of a test particle. 

Let X. n and X out be the homogeneous solutions to Eq.(3.1) ( i.e., with 5 = 0) 
whose boundary conditions axe given by 


X. 


out 


{ 

I 


,—ikr' 


A. e~ ,ur 

in 

B. e~ ikr ‘ 

tn 


+ B oJ hr ' 


for r* — ► — oo ; 
for r* — ► +oo , 
for r* — ► — oo ; 
for r* — > +oo , 


(3.2) 


where k = u — mw + and u> + is the angular velocity of the horizon. Then the 
solution of Eq.(3.1) can be expressed as 


*('•) = w 


X 


out 


(’••) + X (r')J 




(3.3) 


where W is the Wronskian of X. and X 

in out’ 


dX t 

W( u>) = X. out 


tn 


dr* 


-X 


dX. 


out 


dr 


— = 2iwA. (u>) . 

♦ in' ' 


(3.4) 


Since we are interested in the radiation emitted to infinity, we only need the form 
of X at r* — ► +oo, which becomes 


X 


g S(“) e iwr 
2t’u>.4.. (u;) 

in' ' 


G (“0 


oo 

-/■ 


X. Sdr 

tn 




— oo 


(3.5) 


Let h + and h x be amplitudes of the usual plus and cross modes, respectively, 
of gravitational waves at infinity (t.e., the 1/r part of the metric perturbation) 
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(3.6) 


and let h = h , + ih„. Then h is given in terms of G and A. as 
T x w m 


[»1 


h 



G. (w) e - * wu 

>,m v ' 

2iu>c_(u>)A. . (a;) 

O' * tn;i,m v ' 



where S““ is the a = —2 spin-weighted spheroidal function, u = t — r* is the 
retarded time, c Q = + 2) — 12(a 2 u> 2 — mow) — 12iMu and the mode 

indices ( l,m ) tire recovered to avoid ambiguity. 


Let us now evaluate the above integral formula for h. Although it can only 
be done numerically in general, the late time behavior of h can be approximately 
evaluated if one notes the fact that it is dominated by QNMs which are excited by 
a test particle. Furthermore, the dominant contribution comes from those QNMs 
which are least damped. In the present case, they are the ones we obtained in 
the previous section. 


Assuming that the integrant vanishes at large |u;| on the lower half complex 
plane of u; (which should be true for large positive u), the integral in Eq.(3.6) is 
given approximately by the sum of residues of the integrant at poles of 1/A. n (u) 
which are the QNM frequencies; 


fc ~-; E 

l,m=±l *» \ «n v J / 


87rCr f K) Sf n 

c 0 ( w n) w n y/ 2 ir 


e w»« 


(3.7) 


The residues Res (l/A in (u;);a> n ) can be calculated from Eq.(2.4). Inserting the 
result into Eq.(3.7) and approximating all u> n by u> c (= u> c /M) except z n in the 
residues and the exponent of e~ lUnU , we obtain 


h~-- ^ e - ‘^° c *P«(i/2+»(2*«-0) 

r l,m=±l 

47 tG. (uj c ) Sf* 

^ *» m »t m g — IWcti-l -ITTlip 

c o (^c)^c y/2it 
2n7r — 8 


sinh 2n6 


6 cosh 7r(2u> c + 8) 


^2 exp 


26 


(i + *(2w c - «)) - /3uexp 


(3.8) 
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where 


VjJ c M 

Note that the real part of (3 is positive definite (see Eq.(2.9) ). 

Now we evaluate the sum over n in the above equation. Let N be a large 
positive integer and consider the sum over n > N. An inspection of Eq.(3.8) 
shows that this is equivalent to replacing 6 Q by d Q — 2Nit and summing over all 
n > 0. Furthermore, since the important contribution comes from large n, the 
sum over n can be replaced by an integral; 


T Til T / 1 \ nir 0o— 2 Nx 

^2 exp [ — J \ 2 + ^ 2<2>c ~ S y ~ e X e 26 P u 

00 

6 t 8 

— ► — / dx exp [—«® — uue~ x ] = —F(k, v\ u ) , 

TT J IT 


(3.9) 


where 


* = j + *( 2< ^e ~ 0 » 

6n-2Nx 

v = fle 6 . 


The evaluation of the function F(k, u ) is given in Appendix B. It is expressed 
in terms of the incomplete gamma function as 


F(k,v\ u) = ( vu ) * 7 (ac, uu) , 


(3.10) 


where note that 7 (/c, oo) = T(ac). Using Eq.(3.10), we finally arrive at the wave 
form at large u\ 






cosh ir(2u) c + 6) 


.. M\ 4g,„(u, c ) Sfc 


(.. M 


(3.11) 


WcU+*m^ 


, c 0 K)£ c sfa 

where the phase of zu> c should be taken as 7r/2 if u c > 0 and —ir/2 if u> c < 0. 
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The most interesting feature of the wave form is that |h| <x u ~ 1 ^ 2 for u — ♦ oo. 
This implies that the integrated flux of gravitational radiation diverges logarith- 
mically for u — * oo. Thus the extreme Kerr black hole seems marginally unstable 
as suggested by Detweiler! 111 However, the above analysis neglects the presence 
of a test particle and the effect of radiation reaction. Although there is no estab- 
lished method for computing the radiation reaction, one may estimate the effect 
by appealing to the energy and angular momentum conservation laws. 


Let (M, M 2 ) be the mass and angular momentum, respectively, of the initial 
(extreme Kerr) black hole, (er, j) be those of a test particle, and (E r , J r ) be those 
of gravitational radiation emitted to infinity. Then the final mass and angular 
momentum of the black hole, (M^, J^) are given by 


Mj — M + e — E r ; 



(3.12) 


where it is assumed that there is no incoming radiation from infinity. Inserting 
the above into the relation J ^ = Mjqj, solving for A q (= q^ — 1), and retaining 


the terms up to 0[(e/M ) 2 ] we obtain 

_ 2 Me — j e 2 2 ME r — J r 

q ~ ~(M + e) 2 ~M* + M 2 


+ 0 



(3.13) 


It is known that any particle which can fall into the extreme Kerr black hole 
must have e > y/(2M)i 141 Hence if e — j/(2M) = 0(c), the first term dominates 
over the rest and one has A q < 0. However, if c — j /(2M) <C c, all the terms 
may contribute equally to A q. To consider such a case, let us parametrize j as 


j = 2 Me (l - i±) . 
where b is a positive constant of order unity. Then 


A „ c 2 2 ME r - J. 
Ag- (1+24) m! + 


(3.14) 


(3.15) 


10 



It is apparent that the contribution of the test particle is still negative. On the 
other hand, it is not at all clear if the second term is also negative definite, or if it 
is always smaller in magnitude than the first term. To investigate this problem, 
note that E r and J r themselves cannot be evaluated individually, since they are 
logarithmically divergent. Instead, we evaluate the energy flux dE r /dii and the 
angular momentum flux dJ r /du. We have 


du M 2 \ du du ) 


(3.16) 


where dq r /du means the part of dq/du which is due to radiation. 

The formulas for the radiated energy flux and angular momentum flux are 




i_ [ ,L 
du 16tt J 1 

r— ► oo 

*— IS /*■(*») '•"- 


(3.17) 


r— *>oo 


where h = 8h/du and h 1 = dh/d<p. Though the precise evaluation of dE r /du 
and dJ r / du can be done only if the motion of a test particle is specified and only 
by a careful numerical analysis, a qualitative estimate can be done by using the 
result obtained in §3, assuming that the QNMs give the dominant contribution 
to radiation. From Eq.(3.11), for each m = ±1 component, we find 


h = —iu) c h ( 1 + t 
h' = imh , 


(■ 


7 (/c + 1, i/u) 


tu; c u 7 (/e, uu) 


)■ 


(3.18) 


where we have suppressed the indices ( l,m ) for h. Inserting the above into 
Eqs.(3.17), we obtain 
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where only the terms that are important at large u have been retained. Then 
noting that u> c = u> c /M — /(2 M), Eq.(3.16) gives 



m{m — 8) 
16 ttM 2 u 



r—> oo 


(3.20) 


Now let us evaluate the sign of dq T /du. To do so, it should be reminded that 
the Teukolsky equation has the symmetry such that 


<J(u>,l,ro) = Q(-w,l,-m), 


(3.21) 


for any physical amplitude. In particular, this implies 




where G { m is the amplitude appearing in Eq.(3.11). Hence |h Jm | 2 and \h { _ m | 2 
are not independent and one has to sum the both terms if one wants to know the 
sign of {dq T /du) l . From Eq.(3.11), we find 



16 sinh 2 (27r£) e T (^) 
l cosh 7 t(1 ^ 8) cosh 2 7 r(l ± 8) 



(3.22) 


where the upper and lower signs correspond to m = l and m = —l, respectively. 
Hence ( dq r /du) l is given by 


\duJi \du) lt \du) ,_i 

M / e \ 2 I ^|,| I sinh 2 (27r£) e~ K ( l ~ s ) 

u 2 \m) c q tt cosh7r(I — 8) cosh 2 7r(I + 8) 


x 


(l-6) + {l + 8) 


cosh7r(I + 8) 
cosh7r(/ — 8) 



(3.23) 


Thus the sign of ( dq r /du) l depends on the magnitude of 6. 


12 



I 


For / < 6, one can evaluate 8 approximately by using Table 1 of Press and 
Teukolsky. 1 * 71 One finds S l=2 « 2.05 > 2, while 6 < / for / > 3. The latter 
fact implies ( dq T /du )j is positive for / > 3. As for l = 2, a direct evaluation of 
the square bracket term in the last line of Eq.(3.23) shows it is also positive. 
Hence we conclude that dq r /du is positive, i.e., the radiation reaction works in 
favor of the increase in q. We mention that since dq r / du <x u~ 2 , its integral 
A q r = J°° (dq r / du) du is finite, hence can be computed in principle if the test 
particle trajectory is given. Then the result should take the form, 

A - £2 

A?r ~~ n M 2 ’ 

where rj is the efficiency factor which is a function of parameters of the trajectory. 
For a < 1 it has been shown that A q r is negative for all infalling orbits! 41 Thus 
the above result shows a completely distinct nature of an extreme Kerr black 
hole. 



4. Conclusions 

In the lowest order of e/M, it is known that a test particle falling into an 
extreme Kerr black hole cannot increase the value of a/M (= q ) above unity! 1#1 no 
matter what the energy e and the angular momentum j of the test particle are, 
provided e/M 1 and j/M 2 <£. 1. This fact is in accordance with the cosmic 
censorship hypothesis 1141 and with Hawking’s area theorem! 1 ’ 1 since the area of 
the horizon A = SirMr + can increase only if q decreases, if it is unity in the 
beginning. Hence as long as the contribution of the test particle to the change of 
q is O(efM), we have A q < 0. In the case of purely gravitational perturbations 
without a test particle, it is known that the area theorem holds for any value of 
q < l!** 1 Taking the limit q —* 1 from below, one can extend this result and show 
that dq/du < 0 holds for q = 1. However, for certain test particle trajectories 
(j=2Me), the contribution to A q can be 0[{e / M) 2 )]. In such a case, the result 
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in §3 shows that the energy and angular momentum radiated away, which are of 
0[(e/M) 2 )], is important to determine the value of q. If q decreases eventually, 
the black hole will settle down to a stable non-extreme Kerr black hole and the 
logarithmic divergence of the total flux would not have any physical significance. 
On the other hand, if the result is opposite, there arises a possibility of violating 
the cosmic censorship, or at least a possibility of maintaining the value of q very 
close to unity for sufficiently long time. Given the evidences that the cosmic 
censorship does not fail in the both cases of the test particle argument and 
purely gravitational perturbations, perhaps we should be skeptical about the 
former possibility. 

As we noted in the introduction the total angular momentum of the coalescing 
binary neutron stars before the collapse is much larger than that of the extreme 
Kerr black hole. In a sense one can say that the total angular momentum of 
the extreme Kerr black hole is the lowest angular momentum for the binary. 
Unless the system loses its angular momentum down to the angular momentum 
of the extreme Kerr black hole, the system cannot settle down to a single Kerr 
black hole. The system before the final Kerr black hole probably consists of a 
Kerr black hole in the center and a radiating envelope around it. So there may 
exist a case in which the value of q can remain sufficiently close to unity during 
the radiating process and consequently the efficiency of gravitational radiation 
can be very high. Pushing this thought further boldly, one may speculate that 
during the gravitational collapse of a rapidly rotating star or a binary system, 
an extreme Kerr black hole is effectively formed and is kept extreme throughout 
the dynamical stage, emitting a substantial fraction of the rest mass energy in 
the form of gravitational waves. 

If the above conjecture is correct, we cam make some astrophysical impli- 
cations of our results. Clark amd Eardley^ 0 ' studied the evolution of a binary 
neutron star system consists of neutron stars of mass 1.3 Mq and 0.8 Mq and 
estimated the gravitational waves emitted in the stable mass stripping phase. 
They obtained the total energy of the gravitational waves of 6 x 10 52 ergs which 
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is 1.5% of the rest mass. Recently Oohara and N akamura 13 1 1 p erformed Newto- 
nian 3D hydrodynamics calculations of coalescence of binary neutron stars of 
each mass 1.4Af Q . They are using 140 3 grids in the (x,y,z) coordinates. They 
began the simulation when the two neutron stars just touch with each other and 
continued the simulation up to t = 2 msec. The total radiated energy estimated 
by using the quadrupole formula was 2.6% of the rest mass in 2 msec. This ef- 
ficiency can be compared with that extrapolated from perturbation calculations 
for a particle falling into a Kerr black hole with a < 0.99 M , which can be as high 
as 9% in an extreme case! 4 ’ 7 ’' 1 As noted in Introduction, a plausible final destiny 
for coalescence of binary neutron stars is an extreme Kerr black hole. So the 
above conjecture suggests that the efficiency of the gravitational waves is much 
larger than the Newtonian simulations due to the contribution from the QNMs 
for the extreme Kerr black hole. So there is a possibility that ^ 10% efficiency 
is achieved in reality. 

Finally to confirm all the theoretical conjectures in this paper, the construc- 
tion of a 3D fully general relativistic code is strongly needed. 
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APPENDIX A 


Here, we derive the relation between the amplitudes A. n and Z. n . For con- 
venience, we first list some of the coordinate variables which are frequently used 
here and in the text: 


dr * = 


2 , 2 
r + a 


dr, 


dv = dt + dr* , 


dip = dip + —dr , 


A = r 2 — 2Mr + a 2 , 
K = (r 2 + o 2 )w — am , 



(A.l) 


where t, r and p are the standard Boyer-Lindquist coordinates! 3 * 1 Although we 
concentrate on the case a — M in this paper, the following discussion is valid 
even for a M as long as ( M — a)/M 1. ' 

The radial function used in Appendix A of Teukolsky and Press 1131 is of fl, 
in the coordinates ( v , <p), which we denote by R t . The radial function of the 
standard Teukolsky equation 1331 is of T # in the coordinates (t, < p), which we denote 
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by R t . The variables fl t and T, are related to each other by 




Hence the corresponding radial functions are related to each other by 


*-.<*(!)**.«-* •'a*- 


Further, for |s| = 2, there is a simple relation between R 2 and iZ_ 2 ; 

B.,A J « . 

From Eqs.(A.3) and (A.4), one has 


(A.) •-*'**. 


(A.2) 


(A.3) 


(A.4) 


(A.5) 


where f Q is a constant which depends on the normalization of R_ 2 and R_ 2 - l n 
what follows, we denote R_ 2 by R and R_ 2 by R for notational simplicity. 

From Eqs.(2.9) and (2.13) of Sasaki and Nakamura 1 ** 1 , our radial function X 
for the regularized equation is related to R by 






(A.6) 


where /, g and h are, to an extent, arbitrary functions of r except for cer- 
tain boundary conditions they must satisfy, but for definiteness we choose / = 
h =const. and g = (r 2 + a 2 )/r 2 as given in Appendix B of Ref.15. Then X is 
expressed in terms of R as 

* = W + (£ (£*) «- 2i / ■ '**) << S '** . (A-?) 

where /j is a constant. 
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The asymptotic behavior of the M in” mode solution for X (denoted by .X. 
in the text, but here we omit the suffix) is 

X -* l . . . . + (A.8) 

1 A. e tur + A . e twr for r — » oo , 

t. in out ’ 


and that for R is 

( 1 for r — ► r + ; 

R ~* { 2 in (r/r + ) s + Z mt ( r + / r ) c2 '“'' for r — * oo , 


(A.9) 


where r + = Af + \/M 2 — a 2 is the radius of the horizon and u> + = a/ (2 Mr + ) is the 
angular frequency of the horizon. Equations (A.8) and (A.9) fix the normalization 
of the radial functions. 


Now, from Eqn.(A.7) ~ (A.9) and using the solution given in Appendix A of 
Ref.13, it is straightforward to show that A- n is related to Z. n by 

_2 

A. =>/8 -rzZ. , (A. 10) 

tn T in ’ ' ' 

which we quoted in Eq.(2.3). 


APPENDIX B 


Here we evaluate the function F defined in Eq.(3.9); 


OO 


F(k, v ; u) = J dx exp [- 


■kx — vue *] 


Expanding the integrant, we obtain 


00 

F(u)= jdxe- KX ^2- 

l n=0 


(~l) w (i/u) w e~ 
n! 




n=0 

=(t/u) _K 7(/c, i/u), 


(B.l) 


(B.2) 


19 


where 7 (/c, z) is the incomplete gamma function. In particular, we have 


F(u) — ► (uu) *r(«) 


(B.3) 
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